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ABSTRACT 


The  formulation  of  a  one-dimensional  analytical  model 
representing  steady  state  operation  of  a  vertically  oriented 
gas  loaded  variable  conductance  heat  pipe  is  presented.  The 
model  includes  not  only  the  effects  of  binary  mass  diffusion 
and  axial  pipe- wall  conduction,  but  also  gravitational 
effects  which  exist  at  the  vapor-gas  interface  region  when 
the  molecular  weights  of  the  working  fluid  and  the  noncon- 
r’ansible  gas  differ  significantly.  Analytical  expressions 
for  the  transport  processes  {mass,  momentum,  and  energy) 
along  with  equations  of  state  result  in  a  system  of  ten 
equations  which  describe  steady  state  operation  of  the  device. 
Combining  these  expressions  yields  a  system  of  two  simul¬ 
taneous  highly  non-linear  differential  equations.  These 
equations  are  non-dimens ionali zed  and  set  up  for  numerical 
integration  using  a  variable  step  fourth  order  Runge  Kutta 
method  as  part  of  IBM's  Continuous  System  Modeling  Program 
(CSMP)  . 
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I.  INTRODUCTION 

A.  BACKGROUND 

The  heat  pipe  is  a  heat  transfer  device  capable  of  trans¬ 
ferring  large  quantities  of  heat  under  nearly  isothermal 
operating  conditions.  This  device  is  basically  a  liquid-vapor 
two  phase  system  which  transfers  heat  by  evaporation  and 
condensation  of  a  working  fluid  at  opposing  ends  of  a  tube. 
There  are  no  moving  parts  in  a  heat  pipe.  To  sustain  opera¬ 
tion  the  condensed  working  fluid  is  transferred  from  the 
condenser  to  the  evaporator  section  of  the  device  by  capillary 
pumping  action  through  a  wick  structure  lining  the  heat  pipe 
inner  wall. 

B.  CONVENTIONAL  HEAT  PIPE  PRINCIPLES 

The  conventional  heat  pipe  in  its  simplest  form  consists 
of  a  hollow  tube  with  wicking  material  lining  its  inner  sur¬ 
face  as  shown  in  Figure  1.  This  device  is  evacuated,  then 
filled  with  an  appropriate  amount  of  a  working  fluid. 

In  operation,  heat  at  a  given  rate  from  a  source  is  ab¬ 
sorbed  at  one  end  of  the  pipe  evaporating  the  fluid  from  the 
wick.  The  vapor  then  flows  along  the  pipe  to  the  condenser 
end  where  heat  is  rejected  to  a  sink  and  the  vapor  is  condensed 
into  the  wick.  Circulation  of  the  working  fluid  is  maintained 
by  a  capillary  head  in  the  wick  which  must  exactly  balance 
pressure  drops  in  the  liquid  and  vapor,  as  well  as  any  body 
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Figure  1.  Conventional  Heat  Pipe 


forces  present.  The  heat  load  alone  controls  the  satura¬ 
tion  temperature  and  pressure  of  the  working  fluid,  which 
become  the  operating  temperature  and  pressure  of  the  pipe. 

Since  the  transfer  of  heat  is  effected  by  adding  and 
removing  the  heat  of  vaporization  at  saturation  conditions, 
the  heat  pipe  is  nearly  isothermal  along  its  length.  How¬ 
ever,  a  change  in  the  heat  load  will  automatically  change 
the  operating  temperature  (and  pressure)  of  the  pipe,  given 
a  constant  sink  temperature.  This  characteristic  is  unde¬ 
sirable  for  many  applications  where  a  constant  source  tem¬ 
perature  is  required. 

C.  VARIABLE  CONDUCTANCE  HEAT  PIPE  PRINCIPLES 

1.  General 

If  a  noncondensible  gas  is  introduced  into  a  heat 
pipe  along  with  the  working  fluid,  the  heat  pipe  becomes  a 
passively  controlled  device  which  can  operate  at  various  heat 
loads  while  maintaining  a  more  nearly  constant  source 
temperature. 

During  operation  the  working  fluid  behaves  as  described 
above;  however,  the  noncondensible  gas  is  entrained  in  the 
vapor  flow  and  forms  a  gas  plug  at  the  condenser  end  of  the 
pipe.  This  plug  acts  as  a  diffusion  barrier  to  the  flowing 
vapor  so  that  the  effective  condenser  length,  and  therefore 
condenser  heat  transfer  area,  is  reduced.  For  a  given  heat 
load  and  non condensible  gas  loading,  the  pipe  will  operate 
at  a  specific  temperature.  If  the  heat  load  is  increased  the 
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saturation  temperature  and  partial  pressure  of  the  working 
fluid,  as  well  as  total  pressure  in  the  pipe,  tend-  to  in¬ 
crease.  In  response,  the  gas  plug  is  compressed  which 
increases  the  effective  condenser  heat  transfer  area,  thus 
reducing  the  system  temperature  rise.  A  thorough  presentation 
on  the  theory  and  design  of  variable  conductance  heat  pipes 
has  been  provided  by  Marcus  [11. 

2.  Flat  Front  Theory 

Bienert  [2]  developed  the  earliest  model  used  to 
examine  the  operating  characteristics  of  a  gas  loaded  heat 
pipe.  This  flat  front  theory  assumed  that  the  noncondensible 
gas  was  in  the  form  of  a  plug  with  a  sharply  defined  vapor- 
gas  interface  perpendicular  to  the  axis  of  the  pipe.  This 
simple  model  neglected  axial  conduction  in  the  heat  pipe  wall 
and  the  diffusion  of  vapor  into  the  stagnant  gas  plug.  This 
early  model  provided  much  insight  into  the  operating  charac¬ 
teristics  of  gas  loaded  heat  pipes.  However,  experimental 
results  [1]  have  shown  the  model  to  be  inaccurate,  primarily 
because  it  neglects  axial  conduction  and  diffusion. 

3.  Diffuse  Front  Theory 

A  more  accurate  model  of  the  gas  loaded  heat  pipe  was 
provided  by  Edwards  and  Marcus  [3].  The  model,  referred  to 
as  the  'diffuse  front  theory',  includes  effects  of  binary 
diffusion  between  the  working  fluid  and  the  noncondensible  gas 
as  well  as  axial  conduction  in  the  pipe  wall.  Marcus  reported 
[1]  good  correlation  between  experimental  results  and  a  com¬ 
puter  program  based  upon  the  diffuse  front  theory. 


4.  Gravitational  Effects 

Neither  the  flat  front  theory  nor  the  diffuse  front 
theory  included  body  force  effects  on  the  vapor-gas  inter¬ 
face  region.  This  simplification  was  probably  made  because 
much  of  the  early  work  with  heat  pipes  was  related  to  space 
program  applications. 

Kelleher  and  Batts  [4]  reported  convincing  experimen¬ 
tal  results  indicating  that  gravity  effects  can  significantly 
distort  the  axial  temperature  profile  when  the  working  fluid 
and  noncondensible  gas  are  of  significantly  different  molecular 
weights. 

D.  OBJECTIVE 

The  objective  of  this  research  was  to  develop  a  one-dimen¬ 
sional  analytical  model  describing  steady  state  operation  of 
a  gas  loaded  variable  conductance  heat  pipe.  The  model  was 
to  include  not  only  the  effects  of  axial  conduction  and  binary 
mass  diffusion,  but  also  gravitational  effects  on  the  vapor- 
gas  interface  region.  An  additional  purpose  of  this  effort 
was  to  transform  the  analytical  equations  into  a  set  of  non¬ 
dimens  ionalized  simultaneous  differential  equations  which 
could  readily  be  solved  numerically  using  IBM's  Continuous 
System  Modeling  Program  III  [5]. 
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II.  FORMULATION  OF  ANALYTIC  EQUATIONS 

In  formulating  the  analytic  equations  describing  steady 
state  heat  pipe  operation  the  differential  control  volume 
in  Figure  2  was  developed  to  represent  a  section  of  the  heat 
pipe  condenser.  As  implied  in  Figure  2,  only  vertically 
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Figure  2.  Differential  Volume  Element 


oriented  heat  pipes  are  considered  so  that  the  transport 
processes  can  be  treated  as  one  dimensional. 


A.  SIMPLIFYING  ASSUMPTIONS 

The  following  assumptions  were  made  in  formulating  an 
analytical  model  using  the  various  transport  processes: 

1.  The  heat,  mass,  and  momentum  transport  processes  were 
assumed  to  be  one-dimensional  in  the  axial  direction  within 
the  vapor  space.  This  assumption  restricted  this  analytical 
model  to  the  specific  case  of  a  vertically  oriented  heat  pipe. 
In  any  other  orientation,  gravitational  effects  should  be 
expected  rc  cause  highly  three  dimensional  spatial  variations 
when  the  molecular  weights  of  the  working  fluid  and  noncon¬ 
densible  gas  differ  substantially .  This  effect  was  demon¬ 
strated  experimentally  by  Kelleher  and  Batts  [4], 

2.  The  working  fluid  and  the  noncondensible  gas  were  treated 
as  thermodynamically  independent  substances. 

3.  Steady  state  conditions  were  assumed. 

4.  The  working  fluid  was  considered  to  be  at  saturation 
conditions  at  all  times . 

5.  Thermal  resistance  of  the  wick-wall  structure,  in  the 
radial  direction,  was  considered  to  be  negligible  compared 
to  the  thermal  resistance  outside  of  the  pipe. 

6.  The  noncondensible  gas  was  treated  as  a  perfect  gas. 


15  “  pa  +  pb 


(II. 5) 


since  species  'a'  and  species  'b'  occupy  the  same  volume  as 
independent  substances. 

Then  equation  II. 2  becomes. 


Aa  £  +  fe(Vv)  -  fe(pDab  -ar>  ■  0 


(II. 6) 


which  represents  the  conservation  of  mass  for  species  'a'. 

2 .  Species  'b*  (noncondensible  gas) 

A  similar  equation  was  written  for  species  'b'  noting 
that  since  species  'b*  is  a  noncondensible  gas,  m^  e  0. 

Then, 


( w,  pv)  -  i-(pD.  =  0 

dxv  bM  '  dxVM  ba  dx/ 


(II. 7) 


where  the  mass  fraction  of  species  'b'  is  defined  as. 


“b  =  pb/p 


(II. 8) 


3.  Mass  Continuity  of  the  Mixture 

Adding  the  equations  for  species  ’a'  and  species  * b ' 
(III. 6  and  III. 7)  gave  the  continuity  equation  for  the 
mixture , 


_  j  ,  du  die. 


(II. 9) 


where  D  =  is  used  for  simplicity  of  notation. 
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Recalling  equations  II. 4,  II. 5,  and  II. 8, 


(0  +  03.  =  1 

ci  D 


which  reduced  equation  I I. 9  for  mass  continuity  of  the  mix¬ 
ture  to. 


m  ?^  +  -5—  (pv)  =  0 

a  Ac  dx 


(II. 10) 


CONSERVATION  OF  MOMENTUM 


1.  Mixture  of  Species  'a*  and  Species  'b 


Conservation  of  linear  momentum  for  the  mixture,  with 
steady  state  conditions,  was  written  for  the  vapor  space 
differential  element  of  Figure  3  as. 


If  =  /  •  dA) 


where,  dA  is  the  outward  directed  surface  area  element  and, 

— V  «  • 

F  is  the  Siam  of  all  external  forces  acting  on  the  fluid  m 

the  control  volume.  For  the  one-dimensional  case  F  =  F  . 

x 

The  sum  of  the  body  forces  and  surface  forces  was  seen  from 
Figure  3  to  be. 


If  =  p  g  Ac  dx  -  x  Fv  dx  -  dp  Ac  (II. 11) 


and  the  surface  integral,  representing  the  net  efflux  of 
linear  momentum  in  the  axial  direction  from  the  control 
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volume,  was  expanded  to  give: 


i3£Ss 


CS  ,  X 


/  =  -  pv  Ac  +  ( v  +  dv)  ( pv  +  d(  pv)  )  Ac  +  v  ft  Pvdx 

9  9 


=  Ac  pvdv  +  Acvd  {  p v)  +  v  ih  Pvdx 

9  9 


(11.12) 


equaling  equations  II. 11  and  11.12,  and  rearranging  yielded. 


pg 


.  x  Pv  -  |E  =  nV 
Ac  dx 


(11.13) 


For  steady  state  conditions  the  noncondensible  gas  remains 
stationary  with  respect  to  the  coordinates  fixed  to  the  heat 
pipe.  Therefore  the  net  flux  of  species'b'  is  identically 
zero  (i.e.,  pbvb  =  0  since  v^  =  0) .  Noting  that 

cv  =  pava+pbvb'  “d  vb  5  °< 


pv  =  pv 

d  d 


or 


P  v 

V  =  —  v  =  - 

a  p  u) 
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Replacing  v  in  equation  11,13  and  rearranging  resulted  in 
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the  momentum  equation  for  the  mixture. 


d  ,  2,  Pv  .  '“a  Pv  , 

dx  dx  pV  Ac  u  Ac 
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(11.14) 
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where  the  viscous  shear  force  t  is  some  function  of  velocity. 
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For  example,  r  =  ^  v,  if  Poiseuille  flow  is  assumed. 


Species  'b 


Since  the  velocity  of  species  'b'  is  zero,  the  con¬ 


servation  of  momentum  equation  for  species  *b'  could  be 


written  directly  by  inspection  as 


pb  g 


(11.15) 


A  separate  momentum  equation  for  species  'a',  similar 


in  form  to  equation  11.14  could  be  written.  However  the  addi¬ 


tional  equation  would  not  be  useful  since  it  would  not  be 


linearly  independent  of  the  two  momentum  equations  already 


presented. 


D.  ENERGY  BALANCE  FOR  PIPE  WALL 


In  developing  a  heat  pipe  energy  balance  the  wick  and 


pipe  wall  were  considered  together  as  a  composite.  Recall 


that  thermal  resistance  of  the  wick-wall  structure,  in  the 


radial  direction,  was  assumed  to  be  negligible.  Figure  4 


then  represents  an  appropriate  differential  element  upon 


which  to  develop  the  energy  balance. 


By  the  conservation  of  energy  (referring  to  Figure  4) , 


Q.  =  Q  . 
m  out 


then. 
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Qc  +  m  h-  Pv  dx  =  Qc  +  dQc  +  h  Ps  (T  -  T  )  dx 

3.  £  y  00 


(11.16) 


but,  by  Fourier's  law  of  conduction, 


Qc  = 


"  kA*f 


(11.17) 


Incorporating  equation  11.17  into  equation  11.16  and 
rearranging  yields  the  energy  balance  equation. 


^  .  hPs_(  )  +  hfg  Fv 
.2  k  Ap  '  °°;  k  Ap 


=  0 


(11.18) 


where. 


(a)  the  coefficient  of  thermal  conductivity  (k) 
must  be  selected  to  represent  the  composite 
wick-pipe  wall  structure,  and 

(b)  the  heat  transfer  coefficient,  h,  is  some 
function  of  temperature.  The  correlation 
selected  for  h(T)  depends  upon  the  external 
cooling  conditions. 


E.  EQUATIONS  OF  STATE 
1.  Species  'a' 


Since  saturation  conditions  were  assumed  to  exist 
in  the  working  fluid  at  all  times,  tabulated  saturation 
properties  of  the  working  fluid  can  be  used  to  generate 
functions  for  partial  pressure,  partial  density  and  heat  of 
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vaporization  of  the  working  fluid  as  functions  of  local 
temperature .  Then , 


p  =  Fa  (T) 

CL 


(11.19) 


Pa  =  Sa (T) 
a 


(11.20) 


The  functions  Fa(T) ,  Sa(T)  ,  as  well  as  hf^(T)  repre¬ 
sent  actual  tabulated  saturation  properties  of  the  working 
fluid  which  can  be  readily  used  in  a  computerized  interpo¬ 
lation  scheme.  Alternatively  Fa(T)  ,  Sa(T),  and  h^(T)  could 
be  approximated  by  polynomials  or  other  anlaytic  functions 
selected  to  fit  the  working  fluid  saturation  curves. 

2 .  Species  'b* 

As  a  perfect  gas. 


pb  =  Pbir;T 


(11.21) 


3.  Mixture  Properties 


For  independent  substances. 


=  pa  +  pb 


(11.22) 


and. 


(11.23) 
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Although  equation  11.22  is  valid  based  upon  Dalton's  law 
of  partial  pressures,  it  is  included  here  for  completeness 
only  since  it  is  not  linearly  independent  of  equations 
11.14  and  11.15  developed  from  momentum  considerations. 

F.  SUMMARY  OF  EQUATIONS 

The  following  set  of  ten  independent  equations  in  ten 
unknowns  describes  steady  state  operation  of  a  vertically 
oriented  gas  loaded  heat  pipe. 
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(II. 10) 
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(II. 4) 


(11.23) 


p  =  Fa  (T) 


(11.19) 


P_  =  Sa  (T) 


(11.20) 


The  ten  dependent  variables  v,  p,  p&,  p^/  p»  Pa/  P^/ 

co  ,  T)  are  all  functions  of  axial  position. 
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III.  TRANSFORMATION  TO  COMPUTER  MODEL 


A.  ELIMINATION  OF  VARIABLES 

The  ultimate  goal  was  to  obtain  two  nondimensionalized 
differential  equations  with  nondimensional  temperature  and 
mass  fraction  of  species  'a'  as  the  sole  dependent  variables. 

By  reducing  the  number  of  equations  and  dependent  variables 
analytically  the  number  of  computations  required  for  each 
numerical  integration  step  is  significantly  reduced.  This 
approach  was  taken  to  enhance  the  numerical  accuracy  and 
stability  of  computer  results  for  the  model.  With  this 
objective  in  mind,  the  set  of  ten  equations  in  tan  unknowns 
summarized  in  Section  II. F  was  reduced  to  two  equations  in 
two  unknowns  by  the  following  steps. 

1.  Total  pressure  (p)  appeared  in  equation  11.14  only  and 
therefore  could  not  be  eliminated  as  a  dependent  variable  by 
back  substitution.  Since  total  pressure  as  a  function  of 
axial  position  was  not  of  particular  interest,  p  was  eliminated 
by  assuming  that  the  axial  pressure  gradient  of  total  pressure 
was  negligible  (dp/dx  =  0) .  Since  the  axial  pressure  gradient 
of  total  pressure  is  primarily  created  by  the  hydrostatic 
head  of  the  vapor-gas  mixture,  this  common  assumption  is 
reasonable.  This  assumption  means  that  total  pressure  is 
assumed  to  remain  essentially  constant  throughout  the  vapor- 
gas  space  of  the  pipe  for  a  given  operating  power,  gas  load, 
operating  temperature  and  sink  temperature . 


g>; 


It  is  significant  to  note  that  although  the  axial  gradi¬ 
ent  of  total  pressure  is  negligible,  axial  variations  in 
the  partial  pressures  of  both  working  fluid  vapor  and  noncon¬ 
densible  gas  are  significant.  For  example,  using  methanol 
as  the  working  fluid  at  typical  operating  conditions  results 
in  an  axial  partial  pressure  variation  of  several  atmospheres. 

This  assumption  eliminated  equation  11.14  and  the  variable 

?• 

2.  The  following  simple  substitutions  were  then  sequen¬ 
tially  made  in  the  remaining  equations. 

a.  Equation  11.20  was  used  to  eliminate  p  . 

cl 


Pa  =  Sa(T) 

a 


b.  Equation  11.19  was  used  to  eliminate  p. 


pa  =  Fa (T) 


c.  Equation  11.21  was  used  to  eliminate  p^. 


pb  pb  Mfa  T 


d.  Equation  11.23  was  used  to  eliminate  p^. 


pb  =  p  “  pa 


e.  Equation  II. 4  was  used  to  eliminate  p 


Sa  (T) 


i  x 


3.  The  surviving  four  equations  in  four  unknowns  (ih 

v,T,u)  )  were: 
a 


m 


Pv 


a  Ac 


+ 


d 

dx 


(Sa  v) 


d_,Sa_D_  ^a, 

dx '  co  dx 
a 


0 


Pv  d  ,Sa  v 
Ac  dx '  to  ' 


0 


—  fSa(— 1)  — ] 

dxlbata)a  i;MbJ 


Sa  (—  -  1)  g 
0)  3 
a 


Ps  h 
k  Ap 


(T  -  T  ) 
00 


+ 


h_  Pv 

_lS _ 

k  Ap  'a 


0 


4.  Equation  III. 2  was  then  used  to  eliminate  ma  in 
equations  III.l  and  III. 4. 


ih 


Ac  d  ,Sa  v. 
Pv  dx^  u  ' 


Equation  III.l  then  became: 


d  ,Sa  v,  ,  d  t  d  ,Sa  D  da3a, 

-  as'-i— >  +  s(Sa  v)  -  e(-s—  TET> 

3.  3 


=  o 


which  was  integrated  directly  to  yield. 


sa  vil  -  i)  -  *2-2  =  C 

w  _  u)  dx 

«  3 


where  C  is  a  constant  of  integration.  However,  at  x  = 


(111.1) 

(111. 2) 

(111. 3) 

(111. 4) 
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v  =  0 


and  Sa/u)  =  p 
a  o 


then 


C  =  -  D  p  — 3— — 

00  dx 


(III.  5) 


where  x  =  0  was  assumed  to  be  at  the  condenser- reservoir  end 

dw  ! 


of  the  pipe.  However,  -~  =  0  is  a  boundary  condition  as 


will  be  discussed  later. 


Therefore  C  =  0 


The  remaining  system  of  three  equations  was  then. 


Sa  v(l  -  — )  -  ^JE->  a  =  o 
(d  oj_  dx 


(III. 6) 


^[sa<i4'  X)^TI  =  Sal^'  1)9 


(III. 3) 


d2T  _  Ps  h.  ,  hfg  Ac  d  ,Sa  .  „ 

^2  k  Ap  00  k  Ap  dx'  wa 


(III. 7) 


5.  Equation  III. 6  was  finally  solved  for  velocity  and  used 


to  eliminate  v  in  equation  III. 7  to  give  a  system  of  two 


equations  in  T  and  ooa. 


~[Sa(—  -  l)-p-  T]  =  Sa(—  -  1)  g 
dx  u)  M.  u> 

a  b  a 


(III. 3) 


d2T  Psh,  T)  .  hfg  Ac  d  ,  D  Sa  ^a,  =  Q 

,2  k  Ap  k  Ap  dxLu  (to  -1)  dx  J 

dx  3.  3 


(III. 8) 


To  solve  these  two  equations  simultaneously  it  was  advan¬ 
tageous  to  expand  and  differentiate  equation  III. 3  and  to 
expand  and  rearrange  equation  I II. 8  thereby  generating  a 
system  of  two  second  order  differential  equations  in  ua 
and  T. 

After  differentiating  equation  III. 3  a  considerable  amount 
of  algebra  led  to  the  following  form: 

2  d2M 

d  T f  ,1  ,m  dSa  ,  \  a/T  Sa^ 

T~2  ^uT~  1  (T  -dF  +  Sa)}  '  TT-( - 25 

dx  a  dx  u 

a 


—  ^  w  dSa  d  Sa.  ,dTv  2  ,  2  dT  . a_v  dSj*+ga) 

-  (1“)  <2  -g^  +  T  — jMgj}  +  ("2  dx  1FJ  (T  +Sa) 

a  dT  a) 

a 


dT 


-  2 


TSa-^a^A^.l  lwdSa  dTv  Sa 

“3(_dF)  +  ~1(u  (_dF  dxJ  ~~ 

oj  a  u 

a  a 


du 

_ §:} 

2  dx  1 


(III. 9) 


Equation  III. 8  was  manipulated  into  a  similar  convenient 
format  as: 


d2T 


d  u  Ac  h 


iSL 


D  Sa 


dx 


dx 


_ IP 

2  Ap  k  (u  -1)  u 


.}  =  £S_h  } 

1  k  Ap1  «>' 


Ac  h_ 

+  — , — fa 


du 


Ap  k  (u  —  1)  u  dx 

cl  cl 


du 

{-D  Sa(-^-?)  ( 


_  _ _  -f  ..  r.) 

dx  7  xo)  -1  w' 

St.  a 


i  tn  dD  ,  ~  dSa.  dT, 
+  (Sa  dT  +  D  dT)  dx> 


(III. 10) 


Equations  III. 9  and  III. 10  represent  the  final  form  of  the 
desired  equations  in  a  dimensional  format. 


B .  BOUNDARY  CONDITIONS 

To  numerically  integrate  these  second  order  differential 
equations  four  boundary  conditions  must  be  established  at 
x  =  0.  To  avoid  the  obvious  singularity  when  w  =  1.0  (at 

cl 

the  evaporator  end  of  the  condenser)  integration  must  proceed 
from  the  condenser-reservoir  end  toward  the  evaporator. 
Therefore  x  =  0  is  chosen  to  be  the  condenser-reservoir  end 
of  the  heat  pipe. 

Appropriate  boundary  conditions  are  then  as  follows. 


1.  T  =  T  at  x  =  0 

CO 


(III. 11) 


at  x  = 


(III. 12) 


at  x  = 


(III. 13) 


4.  id  =  (u  >0  at  x  =  0 
a  a 

o 


(III. 14) 


The  selection  of  an  initial  value  of  w  is  somewhat  arbi- 

a 

trary.  However,  the  initial  value  of  co  selected,  in  conjunc- 

d 

tion  with  the  values  of  T  and  T  used  in  the  model,  indirectly 

Q  oo  - 

fix  the  amount  of  noncondensible  gas  and  the  operating  power 
level  of  the  heat  pipe  being  simulated  for  a  particular 
working  fluid  -  noncondensible  gas  pair. 


C.  NONDIMENS IONALIZATION  OP  EQUATIONS 


To  nondimensionalize  equations  III. 9  and  III. 10  the 
following  dimensionless  variables  were  defined. 
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xs 


x/L 


nondimensional  length 
(independent  variable) 


T  -  T 

_  00 

Y1  =  ^ nondimensional  temperature 

e  °° 

where : 

T  =  T(x) 

Tg  =  prescribed  evaporator  temperature 
=  ambient  temperature 

Y2  =  In  u>  logarithm  of  the  mass  fraction  of 

species  'a' 


For  ease  of  notation  and  computer  programming  the  first 
and  second  derivatives  of  the  two  dependent  variaoles  were 
represented  as  follows. 
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The  following  nondimensionalized  functions  of  temperature 
were  then  developed  for  direct  substitution  into  equations 
III. 9  and  III. 10. 


D(Te) 


nondimensional  binary  mass 
diffusion  coefficient 


d  (D)  (Te“T<») 

dT  D (T  ) 
e 


nondimensional  temperature 
gradient  of  the  binary  mass 
diffusion  coefficient 


h (T_J 
fg  e 


nondimensional  heat  of 
vaporization 


Sa(Te) 


d(Sa)  ^Te 
dT  Sa(T  ) 


d2(Sa)  (Te-T-> 

Sa<V 


nondimensional  working 
fluid  density 


nondimensional  temperature 
gradient  of  working  fluid 
density 


nondimensional  curvature  of 
working  fluid  saturation 
density  curve 


G7  = 


h(Te) 


nondimensional  convection 
heat  transfer  coefficient 


Using  the  nondimensional  variables  and  functions  of 
temperature  defined  above,  equations  III. 9  and  III. 10  were 
transformed  into  the  following  nondimensionalized  system  of 
two  equations  in  two  dimensionless  dependent  variables  (Y1 
and  Y2) ,  and  their  first  and  second  derivatives  with  respect 
to  axial  (dimensionless)  length. 


DY3{ (1  -  Exp (Y2) )  ( (K1Y1+1) G5  +  K1G4) }  -  DY4 ( (KlYl+1) G4) 


=  ( 1-Exp (Y2) ) {Y3 (Y3 (2K1G5+ (KlYl+1) G6) +K2G5) } 


+  Y4{2Y3 ( (KlYl+1) G5+K1G4)  -  G4 ( (KlYl+1) Y4-K2) }  (III. 15) 


DY3 (Exp (Y2) -1)  -  DY4 (K3G1G3G4)  =  (Exp ( Y2) -1) K4G7Y1 


+  GlG3Y4(K3(G4(Sga  -  +  «Y3) ) 


where  Kl  through  K4  are  dimensionless  constants: 


(III. 16) 


Kl  =  ( 


K2  =  ( 


K3  =  (- 


T  -  T 
e 


“b  ? L  . 

RT  g  1 
00 

Ac  hfq(Te)  D(Te)  Sa(Te) 
k  Ap  (Te  -  TJ 


(III. 17) 


K4  =  ( 


PS  h(T  )  1/ 


k  Ap 


In  terms  of  the  nondimens ional  variables,  the  boundary 
conditions  (equations  III. 11  through  III. 14)  became: 


1.  Yl  =  0 


2.  Y2  =  Y2q  <  0 


at  XS  =  0 


at  XS  =  0 
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j  a 
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3.  Y3  =  DYl  =  0  at  XS  =  0  (III. 20) 

4.  y4  =  DY2  =  0  at  XS  =  0  (III. 21) 

where  y2  =  In  ui  .  As  noted  previously  with  the  dimensional 

cl 

form  of  this  boundary  condition  (equation  III. 14) ,  the  initial 
value  of  a)  selected  is  arbitrary.  However  the  initial  value 

Si 

of  u)a,  in  conjunction  with  the  specified  values  of  Tg  and  T^, 
indirectly  specifies  the  quantity  of  noncondensible  gas  and 
operating  power  level  of  the  heat  pipe.  Since  0  <  <  1.0, 

the  corresponding  value  of  Y2  must  be  negative. 
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IV.  DISCUSSION  AND  RECOMMENDATIONS 

A.  DISCUSSION 

The  Continuous  System  Modeling  Program  (CSMP) ,  developed 
by  the  IBM  Company  [5] ,  is  a  straightforward  computer  package 
available  to  numerically  solve  equations  III. 15  and  111,16. 

The  CSMP  package  includes  a  number  of  optional  numerical 
integration  methods  including  a  variable  step  fourth-order 
Runge-Kutta  method. 

Although  equations  III. 15  and  III. 16  are  highly  nonlinear 
in  the  dependent  variables  Y1  and  Y2,  this  system  of  equations 
can  be  considered  linear  in  the  second  derivatives  of  Y1  and 
Y2.  Therefore  these  equations  can  be  written  in  the  following 
form: 

(C0EF1) (DY3)  +  (C0EF2) (DY4)  =  Fl  (IV. 1) 

(C0EF3) (DY3)  +  (C0EF4) (DY4)  =  F2  (IV. 2) 

where, 

C0EF1  =  { (1-Exp (Y2) ){ (KlY 1+1) G5+K1G4) }  (IV. 3) 


C0EF2 


- ( (K1Y1+1) G4) 


(IV. 4) 


FI  =  (1-Exp (Y2) ) {Y3 (Y3 (2K1G5+ (KlYl+1) G6) +K2G5) } 


F2  =  (Exp(Y2)-l)K4  G7Y1 +  G1  G3  Y4(K3(G4 


G2Y3  Y4Exp(Y2) 

i.-gj  (Exp(il)-i)  G5Y3)) 


(IV. 8) 


Using  these  definitions,  the  nondimensionalized  equations 
representing  steady  state  heat  pipe  operation  may  be  numerically 
integrated  using  the  variable  step  fourth-order  Runge-Kutta 
option  of  the  CSMP  routine  as  follows. 


Yl  =  /  (DYl)dXS  =  / 


d(Yl) 

dXS 


dXS 


(IV. 9) 


+  Y4{2Y3( (KlYl+1) G5+K1G4) -G4 ( (K1Y1+1) Y4 

-  K2)  } 

f 

«  1 

!  1 

‘  | 

(IV. 5) 

-  -4« 

$ 
I 

COEF3 

=  (Exp  (Y2)  -1) 

(IV. 6) 

I 

f 

a 

1 

1 

1 

C0EF4 

=  -(K3G1G3G4) 

(IV. 7) 

i 

& 

a 


$ 

J 

S 


$ 

tl 


Y2  = 


/  (DY2)  dXS  =  / 


d(Y2) 

dXS 


dXS 


(IV. 10) 


Y3  =  /  (DY3)  dXS  = 


/  dXS 

0  dXS 


(IV.. 11) 


I 

% 


1§ 

1 


•*■  »^v<  »•«•*'*’■- 


1  j2 


Y4  =  /  ( DY 4 )  dXS  =  /  d-  dXS 


(IV. 12) 


0  dXS 


where, 


DYl  =  Y3 


(IV. 13) 


DY2  =  Y4 


(IV. 14) 


and  by  Cramer's  rule, 


DY3  = 


(FI  C0EF4  -  F2  C0EF2) 


( COEFl  COEF4  -  C0EF2  COEF3  ) 


(IV. 15) 


DY4  = 


(F2  COEFl  -  FI  COEF  3) 
(OEF1  COEF 4  -  C0EF2  COEF3  ) 


(IV. 16) 


This  model  can  be  used  to  demonstrate  gravitational  effects 
on  the  vapor-gas  interface  region  by  pro.  ling  axial  profiles 
of  temperature  and  the  mass  concentration  of  the  working 
fluid.  When  the  molecular  weight  of  the  noncondensible  gas 
is  less  than  that  of  the  working  fluid,  the  axial  profiles 
will  show  a  fairly  small  interface  region  between  vapor  and 
gas.  However,  when  the  molecular  weight  of  the  noncondensible 
gas  is  substantially  greater  than  that  of  the  working  fluid, 
the  profiles  will  show  a  distorted  (spread  out)  interface 
region. 

Qualitative  results  can  be  obtained  from  this  model  by 
using  additional  simplifying  assumptions  and  built-in  CSMP 
functions  as  follows: 
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1.  Nominal  constant  values  can  be  selected  to  represent 
the  convection  heat  transfer  coefficient  and  the  binary  mass 
diffusion  coefficient.  For  example. 


h(T) 

=  1.0 

BTU/HR/FT2 

for  free  convection 

D  (T) 

=  0.4 

ft2/hr 

for 

Krypton-Methanol 

D(T) 

=  1.9 

ft2/hr 

for 

Helium- Methanol 

2.  To  generate  values  of  working  fluid  density  (Sa)  and 

heat  of  vaporization  (hf^)  at  any  temperature,  a  CSMP  function 

generator  can  be  used  to  provide  quadratic  interpolation  of 

tabulated  saturation  properties.  Tabulated  thermodynamic 

properties  of  methanol  are  available  in  reference  7.  The 

CSMP  function  "SLOPE"  can  be  used  to  generate  values  of 
2  2 

d(Sa)/dT  and  d  (Sa)/dT  as  functions  of  temperature. 

The  temperature  profile  generated  by  this  model  for  a 
given  vertical  heat  pipe  using  a  particular  working  fluid- 
ncncondensible  gas  pair,  with  given  input  values  of  T  , 

T  ,  and  <d  (0)  ,  represents  heat  pipe  operation  at  a  certain 
power  level  and  total  pressure  with  a  specific  load  of  non¬ 
condensible  gas.  The  power  level,  total  pressure,  and  mass 
of  noncondensible  gas  represented  could  be  computed  as 
follows. 
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(IV. 17) 


Mass  of  Gas  (lb) 
m 


1  D  2  L  /  —(l-o)  )  dXS 

4  v  n'  w  a 
0  a 


Power  Level  (BTU/HR)  =  Pv L  /Ah-  dXS  (IV. IS) 

0  a  fg 


Total  Pressure  (psia)  = 


Sa(T  )  R  T 

Fa<T«.)  ‘nf'H11’1 

cl 


“b 

(IV. 19) 


B .  RECOMMENDATIONS 

Based  upon  preliminary  work  using  CSMP,  numerical  results 
from  this  model  could  be  improved  by  the  following  measures. 

1.  Analytical  expressions  are  recommended  in  lieu  of  the 
special  CSMP  functions  to  represent  working  fluid  saturation 
properties  as  functions  of  temperature.  Expressions  for  the 
first  and  second  derivatives  of  density  as  functions  of  tem¬ 
perature  could  then  be  found  analytically  and  used  directly. 
The  resulting  expressions  should  yield  more  accurate  values, 
particularly  for  the  second  derivative  of  density,  than  are 
achievable  through  multiple  use  of  the  CSMP  "SLOPE"  function 
with  tabulated  saturation  values. 

2.  Since  the  binary  mass  diffusion  coefficient  is  a  fairly 
strong  function  of  temperature,  quantitative  results  in  the 
interface  region  could  be  improved  considerably  by  use  of 

an  empirical  expression  for  D(T)  instead  of  a  nominal  constant 
value.  Reference  8  contains  appropriate  empirical  expressions 
for  binary  mass  diffusion  coefficients. 


3.  Quantitative  results  could  also  be  improved  somewhat 
by  using  an  appropriate  empirical  expression  for  the  convec¬ 
tion  heat  transfer  coefficient  as  a  function  of  temperature 
instead  of  a  nominal  constant  value.  Reference  9  provides 

a  variety  of  expressions  for  h(T)  depending  upon  the  convec¬ 
tive  condition  being  simulated. 

4.  With  the  improvements  noted  above,  this  model  could  be 
used  in  an  iterative  loop  to  obtain  preliminary  quantitative 
design  information  concerning  vertically  oriented  terrestrial 
heat  pipes  when  the  molecular  weights  of  the  operating  sub¬ 
stances  are  substantially  different. 

For  example,  if  the  working  fluid- noncondensible  gas  pair 
is  specified,  along  with  the  required  heat  source  (evaporator) 
temperature  range,  ambient  conditions,  and  the  maximum  safe 
total  pressure  for  structural  integrity  of  the  pipe,  then 
equation  IV.  19  could  be  used  to  determine  the  value  of  id  (0) 

a 

to  use  for  boundary  condition  III. 19.  Equation  IV. 17  could 
then  be  used  to  numerically  compute  the  required  mass  of  non¬ 
condensible  gas  to  load  depending  upon  the  heat  pipe  dimen¬ 
sions.  Equation  IV. 18  would  then  yield  the  maximum  heat 
rejection  rate  obtainable  for  these  conditions.  If  this  heat 
rejection  rate  is  not  adequate  for  the  design  requirements, 
the  designer  could  vary  any  of  a  number  of  parameters,  such 
as  operating  temperature,  pipe  dimensions,  total  pressure 
allowed,  substances  used,  etc.,  until  an  adequate  design  is 
found  by  iteration. 
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